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ABSTRACT 

The exact nature of weak Mgll absorbers (those with JJ7(2796) < 0.3 A) is 
a matter of debate, but most are likely related to areas of local star formation 
or supernovae activity outside of giant galaxies. Using 18 QSO spectra obtained 
with the Ultra-Violet Echelle Spectrograph (UVES) on the Very Large Telescope 
(VLT), we have conducted a survey for weak Mgll absorbers at lA < z < 2.4. 
We searched a redshift path length of Az = 8.51, eliminating regions badly 
contaminated by atmospheric absorption so that the survey is close to 100% 
complete to 147(2796) = 0.02 A. We found a total of 9 weak absorbers, yielding a 
number density of absorbers of dN/dz = 1.06±0.12 for 0.02 < 147(2796) < 0.3 A. 
Narayanan et al. (2005) found dN/dz = 1.00 ±0.20 at 0 < z < 0.3 and Churchill 
et al. (1999) found dN/dz = 1.74 ± 0.10 at 0.4 < z < 1.4. Therefore, the 
population of weak Mgll absorbers appears to peak at z ~ 1. We explore the 
expected evolution of the absorber population subject to a changing extragalactic 
background radiation (EBR) from 2 ; = 0.9 to z = 1.78 (the median redshift of 
our survey), and hnd that the result is higher than the observed value. We point 
out that the peak epoch for weak Mgll absorption at z ^ 1 may coincide with 
the peak epoch of global star formation in the dwarf galaxy environment. 
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1. Introduction 


Quasar absorption line systems (QALS) are unique and powerful tools for studying the 
chemical content, kinematics, ionization state, and overall structure of galaxies as well as 
of the intergalactic medium over 0 < z < 6. Mgll absorption is an important feature in 
many QALS, because it relates to metal forming processes. The Mgll AA2796, 2803 doublet 
is a strong transition and is easily detected in low ionization QALS, making it an extremely 
useful tool for probing galaxies and their environments and for studying the structure of 
the intergalactic medium. Weak Mgll absorbers are dehned to be those absorbers with 
Wr{2796) < 0.3 A, and at least some of them represent a unique population(s) of Mgll ab¬ 
sorbers that appears to be much different from stronger absorbers (Rigby et ah 2002; Nestor 
et ah 2005). 

Unlike strong Mgll absorbers, which are typically associated with luminous galaxies 
(within ~ 38h~^{L/ kpc) (Bergeron & Boisse 1991; Bergeron et al. 1992; Le Brun et 
ah 1993; Steidel et al. 1994, 1997; Steidel 1995), weak Mgll absorbers at 0.4 < z < lA are not 
typically found within 50h~^ kpc of such a galaxy (Rigby et al. 2002, but see Churchill et al. 
(2005) for some exceptions), raising questions about their processes of formation. Through 
both statistical means (Churchill et al. 1999; Rigby et al. 2002) and through investigation 
of individual systems (Churchill et al. 2000), weak Mgll absorbers are found to arise in 
sub-Lyman limit systems (15.8 < logA^(Hl) < 16.8 [cm“^]), comprising at least 25% of 
Lya forest clouds in that column density range (Rigby et al. 2002). Weak Mgll absorbers 
are also found to have high metallicity; at least 10% solar, but in some cases even solar or 
supersolar (Rigby et al. 2002; Charlton et al. 2003; Simcoe et al. 2005). This is a striking 
result because, as noted above, luminous galaxies are rarely found within a 50/i“^ kpc impact 
parameter of the quasar. Furthermore, some absorbers have Fell/Mgll ratios that do not 
allow for a-enhancement solely from Type II supernovae, suggesting that metals are also 
produced “m situ” by Type la supernovae and not ejected by nearby, luminous galaxies. It 
is possible that some weak Mgll absorbers arise in areas of intergalactic star formation, or 
are related to supernovae remnants in dwarf galaxies (Rigby et al. 2002). Their incidence is 
therefore likely to be related to global star formation rates and supernovae activity in these 
environments. 

Milutinovic et al. (2005) have found that a flattened geometry for weak Mgll absorbers is 
most consistent with their number statistics and kinematics, and suggest that such structures 
could be part of the cosmic web. Photoionization modeling suggests that the Mgll absorption 
arises in a region ~ 1 — 100 pcs thick (Charlton et al. 2003; Simcoe et al. 2005). This thin 
region that produces Mgll absorption is most likely surrounded by a lower density region that 
produces high ionization Civ absorption, with both regions centered at the same velocity 



(Rigby et al. 2002). Components offset in velocity by tens of km s ^ also often exist in 
Civ and could represent additional surrounding low density filaments (Milutinovic et al. 
2005). 

Previous surveys for weak Mgll absorbers have been conducted a.t 0 < z < 0.3 by 
Narayanan et al. (2005) and at 0.4 < 2 ; < 1.4 by Churchill et al. (1999). Narayanan et al. 
(2005) used E140M Space Telescope Imaging Spectrograph Hubble Space Telescope archive 
data to search for Sill A1260 and Cll A1335 as tracers of Mgll, and found a number den¬ 
sity of absorbers dN/dz = 1.00 ± 0.20. Churchill et al. (1999) used the High Resolution 
Spectrograph on the Keck I Telescope and observed the Mgll AA2796, 2803 doublet directly, 
Ending dN/dz = 1.74 ± 0.10. These values imply that the total cross section for weak 
Mgll absorption is twice that of strong Mgll absorbers (Steidel & Sargent 1992). 

The Endings of Narayanan et al. (2005) for 0 < z < 0.3 weak Mgll absorbers are 
consistent with cosmological evolution of the z ^ 1 weak Mgll absorber population, based 
upon ACDM cosmology (Churchill et al. 1999). However, when the efiect of the evolution 
of the extragalactic background radiation (EBR) on the absorber population is taken into 
account, dN/dz is expected to be higher at ~ 0 than observed. This is partly because 
the less intense ionizing background (compared to ~ 1) leads to an increase in the low- 
ionization Mgll phase, i.e. absorbers that were below the equivalent width threshold at 
z ^ 1 are detectable at z ^ 0. Also, there is a signiEcant contribution at z ^ 0 from 
Mgll absorption in a lower density phase that gave rise only to higher ionization absorbers at 
z ^ 1 (Narayanan et al. 2005). It appears that the processes that lead to weak Mgll absorber 
formation are less active at 2 : ~ 0 then at z ^ 1. Determining the evolution of dN/dz at 
higher redshifts is an important step toward identifying the formation processes of weak 
Mgll absorbers, and thus served as the motivation for this work. 

In order to push observational measurements of dN/dz for weak Mgll absorbers to 
higher redshifts, we surveyed the spectra of 18 QSOs obtained with the Ultra-Violet Echelle 
Spectrograph (UVES) on the Very Large Telescope (VLT). Such a measurement will shed 
light on the nature of the absorber population, and may indicate which factors afiect the 
evolution of the absorbers. The data and survey method are outlined in § 2. Our survey 
results are summarized in § 3. We also simulated the evolution of the = 0.9 absorber 
population back to 2 ; = 1.78, subject to the changing EBR in order to predict the expected 
dN/dz evolution. These results are presented in § 4 and compared with our survey results. 
In § 5 we summarize our Endings and end with a discussion. 
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2. Data and Survey Method 
2.1. VLT Spectra 

The 18 UVES/VLT QSO spectra used for our survey of weak Mgll absorption were 
obtained from the ESO archive. The spectra typically covered a wavelength range of ~ 
3050 — 10075 A, though there were usually gaps in coverage from ^ 5753 — 5841 A and 
~ 8520 — 8665 A (see Figure 1). The resolution was R ~ 45,000, or ~ 6.7 km s“^. Data 
reduction procedures are described in Kim et al. (2004). Continuum htting followed standard 
procedures using the IRAF"^ SFIT task. A cubic spline function was used, typically of order 
7-9. 

Though the spectra were of extremely high quality, there were some complicating issues. 
As already mentioned, there were some small breaks in the coverage in many of the spectra. 
The dense Lya forest made a. z < 0.7 Mgll AA2796, 2803 doublet survey impractical for some 
quasars, but this was not a concern since the Churchill et al. (1999) survey also covered this 
range. There was also signihcant contamination from telluric absorption features, making a 
survey of affected wavelengths impractical, as well. Our general philosophy was to exclude 
troublesome regions rather than to risk missing Mgll doublets in our survey. 

2.2. Survey Method 

For our survey, we hrst assumed that every line detected above a 5a signihcance level 
was Mgll A2796. It was then determined if the corresponding Mgll A2803 line was present 
at the expected location with at least a 2.5a signihcance level. If this was the case, the 
line prohles were checked for similarity in shape by visual inspection, and for consistency 
with a 2:1 optical depth ratio. In all cases some supporting lines at the redshift of our 
Mgll candidate were found (such as Lyo;, Fell, Aim, and/or Civ). Once the detection of a 
Mgll absorber was conhrmed, we measured the rest frame equivalent width. In order to be 
accepted as a weak absorber, a rest-frame equivalent width of 0.02 < IIA(2796) < 0.3 A was 
required. We note that our survey could have extended to a much lower equivalent width 
limit. The 117(2796) = 0.02 A limit was chosen for comparison to previous surveys. This 
corresponds to log A^(Mgll) = 11.72 [cm“^] for a typical Doppler parameter of 6 = 4 km s“^, 
for an unsaturated line. Fortunately, above this limit interference by telluric absorption did 
not affect our survey in our selected wavelength regions (i.e. telluric features were weaker 

"‘IRAF is distributed by the National Optical Astronomy Observatories which are operated by AURA, 
Inc., under contract to the National Science Foundation. 
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than hh^(2796) = 0.02 A). 

For the purpose of the present survey, we focused on absorbers whose redshift fell in 
the range lA < z < 2.4. Several quasars had coverage at z ^ 2.5, but the path length was 
not signihcant at that 2 :. The large wavelength coverage of these spectra also allowed us to 
search for weak Mgll absorbers to signihcantly lower redshifts in most spectra. Therefore, 
we also considered absorbers with redshifts in the range 0.4 < 2 ; < 1.4 for comparison with 
the results of Churchill et ah (1999). 

The excellent S/N of these spectra allowed our survey to be nearly 100% complete if 
we eliminated wavelength ranges with heavy atmospheric absorption. The telluric features 
in the eliminated regions acted as “effective noise” which made unambiguous detection of a 
0.02 < Wr{2796) < 0.3 A Mgll A2796 line at a 5a signihcance level difficult, if not impossible. 
Figure 1 shows the redshift coverage for the Mgll AA2796, 2803 doublet in each of our spectra. 
Because we excluded certain redshifts from our search and because not every QSO spectrum 
covered the Mgll AA2796, 2803 doublet up to 2 : = 2.4, our average redshift does not coincide 
with the center of our redshift range. Considering these factors, the average redshift is given 
by 

^high 

f zdz 

«= E (1) 

included regions I dz 

^low 


In order to verify that we have nearly 100% completeness in the spectral regions that 
we surveyed, we randomly added two artihcial Mgll doublets to each of our spectra, in the 
redshift range 0.4 < < 2.4. Simulated lines were not placed in regions of the spectra that 

were not searched, i.e. the Lya forest, telluric absorption regions, or coverage gaps. The 
artihcial lines were all assigned the minimum equivalent width for our survey, IT)-(2796) = 
0.02 A. Lines were added, beginning with the observed hux and convolving it with the line 
spread function for the simulated line. This allowed us to retain a realistic representation 
of noise in the region, making artihcial lines indistinguishable from real ones. The spectra 
were searched for these artihcial lines in the same manner that we searched for real lines. 
All the artihcial lines were found even though they were at the minimum equivalent width 
for our survey. Therefore, we are conhdent that we ehectively have 100% completeness in 
the regions we searched. 
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3. Results of the Survey 

Survey completeness as a function of redshift was calculated according to the formalism 
given in Steidel & Sargent (1992) and Lanzetta et al. (1987). The total redshift path covered 
in the 18 QSO spectra over the range 1.4 < < 2.4 is given by 

2.4 

Z{Wr,R) = j g{Wr,z,R)dz (2) 

1.4 

where g{Wr, z, R) is the function that gives the number of lines of sight along which the 
Mgll AA2796,2803 doublet could have been detected at a redshift z and with a rest frame 
equivalent width hlA(2796) at greater than or equal to a 5a significance level. R is the 
expected ratio of equivalent widths between Mgll A2796 and Mgll A2803, which, according to 
atomic physics, is two for unsaturated lines. As described in § 2.2, we find 100% completeness 
for all Mgll doublets with 147(2796) > 0.02 A. Therefore, we have the simplest case in which 
Z(0.02 A, 2) = = 8.51, the full path length covered by our survey. 

For the 1.4 < < 2.4 Mgll doublets found in our survey. Table 1 lists the absorber 

redshifts (calculated from the optical depth mean) and the quasar in whose spectrum they 
were detected. Column 3 lists the velocity width of the systems (as defined in Churchill 
& Vogt (2001)), Column 4 the rest frame equivalent width of Mgll A2796, and Column 5 
the Mgll doublet ratio, DR = 147(2796)/1T7(2803). Their Mgll A2796 and Mgll A2803 
profiles are shown in Figures 2-10. We find 9 systems over a redshift path of = 8.51 
at lA < z < 2.4 in the 18 spectra surveyed. Because our survey is 100% complete to 
147(2796) = 0.02 A, we can simply divide the number of systems by the redshift path to 
obtain a redshift number density dN/dz = 1.06±0.12. Figure 11 shows our result, along with 
the results of Narayanan et al. (2005) and Churchill et al. (1999). We also show our survey 
results for the redshift range 0.4 < < 1.4, separated into the same bins as in Churchill et 

al. (1999), for which we find agreement to within l — 2a. Overall, for 0.4 < 2 ; < 1.4, we found 
29 systems in a redshift path of A;^ = 14.5, yielding dN/dz = 2.01 ± 0.14. Churchill et al. 
(1999) found dN/dz = 1.74 ± 0.10 for a somewhat larger redshift path of A^ = 17.2. The 
solid curve shows the no-evolution expectation for dN/dz in a ACDM universe (0^ = 0.3 
and Oa = 0.7) normalized to dN/dz = 1.74 at z = 0.9. 
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4. Expected Evolution of the Weak Mgll Absorber Population Due to the 

Changing EBR 

The number density of weak Mgll absorbers changes with redshift for a variety of 
reasons. Since we expect weak Mgll absorption to arise in small-scale regions, the process 
of hierarchical structure growth could destroy the absorbers through mergers or collapse. 
As the metallicity of the universe increases, we expect to see an increase in the number of 
metal-rich systems. Also, the cumulative effects of supernovae, causing superbubbles and 
superwinds, will redistribute gas in the universe. The ionization state of weak Mgll absorbers 
is also expected to change as the EBR evolves. The EBR is known to change with redshift 
because of the change in the number density of bright quasars and in the rate of global star 
formation (Haardt & Madau 1996, 2001). This changing EBR will have an effect, even on 
an otherwise static population of absorbers, because of their changing ionization state. We 
explore this effect in the following sections. 


4.1. Expected Change in the Ionization State of Mgll Absorbers 

Models by Rigby et ah (2002) and Charlton et ah (2003) suggest that weak Mgll ab¬ 
sorbers consist of a high density region of Mgll and a lower density Civ phase. Typical 
densities are on the order of 0.1 cm“^ for the region that gives rise to the majority of the 
Mgll absorption and 0.001 cm“^ for the Civ region. We conduct a thought experiment 
to explore how these absorbers evolve. Let us assume that a given absorber is stable over 
time, i.e. there is no gravitational collapse, fragmentation, mergers, evaporation, or pressure 
changes. In other words, assume that the absorber is subject to no changes except for the 
evolving EBR. The total hydrogen column density, log A(Htot), the electron number density, 
rte, and the metallicity, logZ/Z©, are all constant. One may ask how the ionization state of 
the absorber will change as a result of the changing EBR, between = 0.9 and = 1.78. 
At = 1.78 the EBR is ~ 0.4 dex more intense than at z = 0.9, though the shape of the 
spectrum is similar. This increase in the photon number density, leads to a change in the 
ionization parameter, [/ = n^/ue- For a given cloud we can fix log A(Htot) and calculate the 
change in logA(Hl) and logA(Mgll). At higher redshift, we find that the increased EBR 
causes the equivalent width to be larger for the high ionization state Civ and smaller for 
the low ionization state Mgll, compared to intermediate redshifts. 
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4.2. Simulated Evolution of Three Systems 

To illustrate the evolution of weak Mgll absorbers, we simulate the appearance at 
2 : = 1.78 of three systems modeled by Charlton et ah (2003) at {z) = 0.9 subject to the 
changing EBR and cosmological evolution. We applied photoionization models using the code 
Cloudy (version 94.00; Ferland et ah 1998) with an EBR of the form prescribed by Haardt & 
Madau (1996, 2001). We adopted their model including quasars and star forming galaxies, 
with an escape fraction of 10% and including absorption by the intergalactic medium. The 
results are presented in Figures 12 - 14. Each of the three systems has only a single cloud 
detected in Mgll, but their CiV absorption prohles differ in strength, kinematic spread, and 
number of clouds. The 2 : ~ 0.9 (Figure 13) and 2 ; ~ 0.6 (Figure 14) systems have offset 
Civ components in addition to those centered on the Mgll. 

Each system was evolved to z = 1.78 subject to the changing EBR. Once the column 
density of each component was found, we created a synthetic, noise-free spectrum. For the 
~ 0.6 system, the Mgll prohle evolves below the equivalent width cutoff for our survey 
(see Figure 14). The overall trend for all three systems is for weak Mgll absorbers to be less 
prevalent at 2 ; = 1.78 than a.t z = 0.9, while Civ absorption, which traces the high ionization 
state, becomes stronger. 


4.3. Expected dN/dz at 2 : = 1.78 Due to the Changing EBR 

We apply a simple physical model to predict dN/dz at z = 1.78 using the results of 
Churchill et ah (1999) as a referemce point. This model only takes the evolving EBR into 
account but is useful for exploring which factors affect the evolution of dN/dz for weak 
Mgll absorbers. Based on the previous measurement of dN/dz at z = 0.9 by Churchill 
et ah (1999), we now consider what dN/dz would be expected at z = 1.78 subject to the 
changing EBR. We first determine the 1W(2796) of a z = 0.9 absorber whose ionization 
state has increased just enough from 2 : = 0.9 to 2 : = 1.78 that it falls below the detection 
threshold of our survey. We compute the column density of Mgll for an absorber with the 
threshold 1W(2796) = 0.02 A to be logiV(Mgll) = 11.7 [cm ^], assuming a typical Doppler 
parameter, b = A km s“^. Models of weak Mgll absorbers (Charlton et ah 2003; Rigby et al. 
2002) have yielded a range of metallicities, —1.0 < \ogZ/ZQ < 0.0, and a range of densities 
10“^ < Ue ;< 10“^ cm“^. To illustrate our method we first choose the values that yield the 
smallest predicted dN/dz value at z = 1.78, XogZ/Z^ = —1.0 and Ue = 10“^ cm“^, and 
solve for the total column density logiV(H)tot = 18.3 [cm“^] subject to the Haardt-Madau 
EBR at 2 : = 1.78, again using Cloudy. We then compute the 1T).(2796) of an absorber with 
this same total column density subject to the EBR at 2 : = 0.9. We hnd 147.(2796) = 0.041 A. 



This is the lowest equivalent width Mgll absorber at z = 0.9 that would have been observed 
as a weak Mgll absorber at z = 1.78. A similar procedure was followed to find an upper 
limit on MA(2796) for absorbers that exceed the 0.3 A threshold for weak absorption at 
= 0.9 but were weak at z = 1.78. Such an absorber will have a total column density 
logA^(H)tot = 21.4 [cm“^], which will produce a 144(2796) = 0.31 A absorption at z = 0.9. 

Once an upper and lower limit on equivalent width at z = 0.9 have been obtained, we 
may integrate over the equivalent width distribution function at this redshift to estimate the 
number of absorbers per unit redshift that we would expect to observe at z = 1.78. The 
equivalent width distribution follows a power law form: 

n(W)dW = CW-^dW (3) 


where C k, 0.4 and 5 = 1.04 (Churchill et ah 1999). The expected ratio of dN/dz at z = 1.78 
to dN/dz at z = 0.9, due only to the changing EBR, is given by 


{dN/ (i^)i.78 
{dN/dz)o.9 


0.31 

/ n{W)dW 

0.041 _ 

0.3 

/ n{W)dW 
0.02 


0.74 


( 4 ) 


Cosmological evolution of the absorber population must also be taken into account. For a 
model with Ha = 0.7 and Qm = 0.3, cosmological evolution alone would predict a factor of 
1.33 more absorbers at z = 1.78. Thus, when we combine evolution due to the EBR and 
cosmology, we find a predicted dN/dz = 1.71 at z = 1.78. 

If we vary the input cloud properties, such that \ogZ/ZQ = 0.0 and Ue = 10“^ cm“^, 
we find a predicted dN/dz = 2.14 at z = 1.78. In that case, the electron density change 
dominates, and the metallicity change only has a small effect. For this estimate we assumed 
that typical weak Mgll doublets are resolved so that 6 = 4 km s“^ is a typical value. There 
may be some unresolved saturation in some of the lines in our sample, since the donblet 
ratios of weak lines are not equal to 2 (see Table 1). Assuming 6 = 2 km s“^, we see 
no significant change in onr predicted dN/dz. Also, 6 = 2 km s“^ implies T < 5800 K, 
somewhat smaller than the equilibrium temperature Cloudy finds for our assumed density 
and metallicity. An even smaller 6, which we consider nnlikely, would result in somewhat 
smaller dN/dz predictions, but would also require a very cold cloud. Thus we believe that 
6 = 4 km s“^ provides a realistic estimate. 

The predicted dN/dz of 1.71 — 2.14 is significantly higher than the observed valne of 
dN/dz = 1.06±0.12 for weak Mgll absorbers at (z) = 1.78. This implies that the population 
of objects that give rise to weak Mgll absorption are becoming more common at z = 0.9 



than they were at z = 1.78. This could occur either through gradual formation of these 
objects, assuming that they are not destroyed, or through an increase in the importance of 
processes that generate a non-static population. 


5. Summary and Discussion 

We have conducted a survey of 18 QSO spectra in the ESO archive obtained with the 
UVES on the VLT, searching for weak Mgll absorbers with 0.02 < 114(2796) < 0.3 A. In 
order to test our survey completeness we randomly added simulated Mgll AA2796, 2803 tran¬ 
sitions to the spectra. These simulated lines were at the low end of the equivalent width 
limit. Because the spectra have high signal-to-noise, the survey is nearly 100% complete in 
the regions that we searched, having eliminated regions strongly affected by telluric absorp¬ 
tion. We found 9 systems in a redshift path length of = 8.51 over the range 1.4 < < 2.4. 

This yields a number density per unit redshift of dN/dz = 1.06 ± 0.12 at (z) = 1.78. For 
comparison, Narayanan et ah (2005) hnd dN/dz = 1.00 ±0.20 at (z) = 0.15 and Churchill et 
ah (1999) hnd dN/dz = 1.74 ± 0.10 at (z) = 0.9. We conclude that the weak Mgll absorber 
population peaks at 2 : ~ 1, and that the objects that generate this type of absorption are 
less common both at higher and at lower redshifts. 

A decrease in dN/dz from 2 : = 1.78 to z = 0.9 would be expected simply on the basis 
of cosmological evolution. The strength of Mgll absorption coming from a population of 
objects is also expected to change subject to the changing EBR. In order to understand 
the balance of these effects, we calculated the expected dN/dz at z = 1.78 based upon the 
known population at z = 0.9. For the purpose of this estimate, we assumed that the absorber 
population was static or that destruction was balanced by the formation processes between 
these two epochs, i.e. we only accounted for evolution due to the changing EBR and due to 
cosmological evolution. Our prediction of dN/dz = 1.71 — 2.14 is markedly higher than our 
observed value. Therefore, we conclude that one of two general scenarios must apply. First, 
the structures that give rise to weak Mgll absorption could be long-lived, and there could 
be a gradual build-up of these structures over time, such as through hierarchical structure 
formation. Second, the structures could be transient, such that the rate of their formation 
process determines dN/dz. An obvious example of such a process is star formation. 

This second possibility was advocated by Narayanan et al. (2005) who considered the 
evolution of dN/dz from 2 : = 0.9 to = 0. In that case, starting with the observed 2 : = 0.9 
population of weak Mgll absorbers, the evolving EBR along with cosmological evolution 
would predict a signihcantly larger population at z = 0 than is observed. This issue was 
complicated by the incidence of a new type of low redshift weak Mgll absorber descended 
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from absorbers from which only CiV absorption was observed at z = 0.9. Nonetheless, the 
trend for an overall decrease in the weak Mgll absorber population is consistent with the 
decrease in the global star formation rate over that redshift interval. This suggests a possible 
relationship between star formation and weak Mgll absorbers, which is also consistent with 
their high metallicities and (in at least some cases) iron abundances (Rigby et ah 2002). 

We should therefore consider how the global star formation rate relates to the results of 
our (z) = 1.78 survey. The global star formation rate is roughly constant over the redshift 
range 1 < z < 4, and then decreases to lower redshifts (see Gabasch et ah (2004) and 
references therein). If production of weak Mgll absorbers was directly correlated with the 
global star formation rate we would expect a dN/dz at z = 1.78 consistent with a constant 
creation rate, subject to the EBR and to cosmological evolution. Our observation of a smaller 
dN/dz than this expectation would suggest that the global star formation rate should be 
smaller at z = 1.78 than at z = 0.9. This is not the case. However, we need to consider that 
weak Mgll absorbers are preferentially found in certain environments. 

Theoretically, in the context of hierarchical structure formation and empirically on the 
basis of the morphology-density relationship, it is expected that star formation peaks later 
in low mass galaxies (Kauffmann et ah 2004). More directly, Bauer et ah (2005) found a 
shift over the range 0 < < 1.5 in the relative contributions of low and high mass galaxies 

to the global star formation rate, with low mass galaxies contibuting more at low redshift. 
Although measurements of the dwarf galaxy contribution have not yet extended to higher 
redshifts, it seems likely that the star formation rate in dwarfs peaks at a lower redshift 
than that for giants. Given the constant global star formation rate at 1 < < 4 (Gabasch 

et al. 2004), if the giants are contributing relatively more at high redshift, the dwarfs must 
contribute less. Given previous suggestions that weak Mgll absorbers are related to dwarfs, 
it is intriguing that the peak epoch of weak Mgll absorbers may coincide with the peak of 
star formation in the dwarf environment. 

Rigby et al. (2002) estimate that weak Mgll absorbers account for approximately 25%- 
100% of Lya forest clouds with 15.8 < logiV(Hl) < 16.8 [cm]“^ at z ^ 1. At first glance it 
seems surprising that the redshift path density for Lya forest clouds would sharply decrease 
with decreasing until = 1.5 {dN/dz oc (1 + ;2)2-47±o.i8^ (Kim et al. 2002; Weymann et al. 
1998), while we find the weak Mgll absorber population to increase over the same interval. 
However, if the weak Mgll absorbers are closely related to star formation in dwarfs, as we 
have suggested, then there is no reason to expect a direct correlation of their evolution with 
the evolution of the Lya forest, the latter which is explained by cosmological expansion. 
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Table 1. 


Properties oi lA < z < 2A Weak Mgll Systems 


QSO 

^abs 

(km s“^) 

1W(2796) 

(A) 

DR 

HE2347-4342 

1.405362 

5.8 ±0.2 

0.076 ±0.001 

1.84 ±0.04 

Q0002-422 

1.446496 

4.9 ±0.2 

0.042 ±0.001 

1.65 ±0.04 

Q0122-380 

1.450109 

28.4 ±0.8 

0.061 ±0.003 

1.34 ±0.12 

HE2217-2818 

1.555845 

30.9 ±0.1 

0.268 ±0.001 

1.44 ±0.01 

HEOOOl-2340 

1.651462 

4.0 ±0.3 

0.070 ±0.001 

1.51 ±0.03 

HE0151-4326 

1.708494 

4.0 ±0.3 

0.027 ±0.001 

1.96 ±0.10 

HE2347-4342 

1.796237 

5.3 ±0.1 

0.146 ±0.001 

1.19 ±0.01 

Q0453-423 

1.858380 

38.6 ±0.4 

0.254 ±0.002 

1.21 ±0.01 

HE0940-1050 

2.174546 

6.9 ±0.4 

0.028 ±0.001 

1.56 ±0.10 
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Q0453-423 
Q0420-388 
Q0329-385 
Q0122-380 
Q0109-3518 
Q0002-422 
PKS2126-158 
PKS1448-232 
PKS0329-255 
PKS0237-23 
HE2347-4342 
HE2217-2818 
HE1341-1020 
HE1158-1843 
HE1122-1648 
HE0940-1050 
HE0151-4326 
HE0001-2340 


0 .5 1 1.5 2 2.5 

Redshift Coverage of Mgll Doublet 


Fig. 1.— The redshift coverage for the Mgll AA2796, 2803 doublet in each of our 18 QSO 
spectra. Our survey covers a total redshift path length of Az = 8.51 for 1.4 < z < 2.4. 
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Fig. 2.— Absorption profiles, plotted in velocity space, for Mgll AA2796,2803 for 
the 117(2796) = 0.076 ± 0.001 A absorber at z = 1.405362 toward HE2347-4342. 
Civ AA1548,1550 is also detected at this redshift. 


1 


0 


1 


0 


Fig. 3.— The HA(2796) = 0.042 ± 0.001 A absorber at z = 1.446496 toward Q0002-422. 
Civ AA1548,1550 is also detected at this redshift. 
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-200 -100 0 100 200 
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Fig. 4.— The 1F^(2796) = 0.061 ± 0.003 A absorber at z = 1.450109 toward Q0122-380. 
SilV AA1393,1402 and CiV AA1548,1550 are also detected at this redshift. 
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Fig. 5.— The 1F^(2796) = 0.268 ± 0.001 A absorber at z = 1.555845 toward HE2217-2818. 
Lya, Sin A1260, Cii A1335, Siiv AA1393,1402, Civ AA1548,1550, Alii A1671, Fell A2344, 
Fell A2383, and Fell A2600 are also detected at this redshift. 
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Fig. 6.— The 1F[.(2796) = 0.070 ± 0.001 A absorber at z = 1.651462 toward HEOOOl- 
2340. Lya, Sim A1207, Sill A1260, Cii A1335, Siiv AA1393,1402, Sin A1527, and 
Civ a a 1548,1550, are also detected at this redshift. 
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Fig. 7.— The 1F[.(2796) = 0.027 ± 0.001 A absorber at z = 1.708494 toward HE0151- 
4326. Lya, Siiv AA1393,1402, and Aim AA1855,1863 are also detected at this redshift. 
The expected locations of many other common features fell in the Lya-forest, making their 
identihcation difficult. The anomaly near the Mgll A2796 line is a residual due to poor sky 
subtraction. 
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Fig. 8.— The 1F[.(2796) = 0.146 ± 0.001 A absorber at z = 1.796237 toward HE2347- 
4342. Lya, Sill A1260, Cii A1335, Siiv AA1393,1402, Sill A1527, Civ AA1548,1550, 
Aim AA1855,1863, Fell A2383, and Fell A2600 are also detected at this redshift. 
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Fig. 9.— The 1F^(2796) = 0.254±0.002 A absorber at z = 1.858380 toward Q0453-423. Lya, 
Sill A1190, Sill A1193, Sim A1207, Sin A1260, Cii A1335, Aim A1671, Aim AA1855,1863, 
Fell A2344, Fell A2374, Fell A2383, and Fell A2600 are also detected at this redshift. 
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Fig. 10.— The 1F^(2796) = 0.028 ± 0.001 A absorber at z = 2.174546 toward HE0940-1050. 
Lya, Sim A1207, Sill A1260, Cll A1335, SilV AA1393,1402, and Civ AA1548,1550 are also 
detected at this redshift. 







22 


N 

T3 


T3 

H—' 

w 

c 

CD 

Q 


CO 

Q_ 


CO 

■Q 

CD 

cr 



Redshift, z 


Fig. 11.— The redshift path density {dN/dz) for our survey over the redshift range 1.4 < 
2 ; < 2.4 is presented as a point with solid error bars. We hnd dN/dz = 1.06 ± 0.12. We also 
present our results for the redshift range 0.4 < 2 ; < 1.4 with solid error bars, for comparison 
to the results of Churchill et al. (1999), that are shown with dotted error bars. The result 
of Narayanan et al. (2005) is shown with dashed error bars. The solid surve shows the 
expected evolution of absorbers based an a ACDM cosmology. The series of black dots 
indicates the range of expected dN/dz for our survey based upon the evolving EBR and a 
ACDM cosmology, assuming a static population. 
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Fig. 12.— The z = 0.8 system towards PG1634+706 (Churchill et al. 1999). The bottom 
panels show the same system evolved to z = 1.78. The left panels show the actual spectra 
superimposed with the full models, and the right panels show the various components from 
the models. This system has a single cloud in Mgll (centered at 0 km s“^) and, although its 
single Civ component arises in a separate phase, it is aligned with the Mgll in velocity. 
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Fig. 13.— The z = 0.9 system towards PG1634+706 (Churchill et al. 1999). The bottom 
panels show the same system evolved to .^ = 1.78. The left panels show the actual spectra 
superimposed with the full models, and the right panels show the various components from 
the models. This system has a single cloud in Mgll (centered at 0 km s“^) but has multiple 
Civ components. 
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Fig. 14.— The 2 ; = 0.6 system towards PG1634+706 (Churchill et al. 1999). The bottom 
panels show the same system evolved to z = 1.78. The left panels show the actual spectra 
superimposed with the full models, and the right panels show the various components from 
the models. This system has a single component in Mgll (centered at 0 km s“^), but has 
multiple components in Civ. The Mgll for this system evolves to be below the detection 
criterion, while the Civ becomes stronger. 


























































































